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Infrared Vortex-State Electrodynamics in Type-II Superconducting Thin Films 
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The vortex-state electrodynamics of i-wave superconductors has been studied by infrared spectroscopy. 
Far-infrared transmission and reflection spectra of superconducting Nbo.5Tio.5N and NbN thin films were 
measured in a magnetic field perpendicular to the film surface, and the optical conductivity was extracted. 
The data show clear reduction of superconducting signature. We consider the vortex state as a two- 
component effective medium of normal cores embedded in a BCS superconductor. The spectral features 
are well explained by the Maxwell-Garnett theory. Our analysis supports the presence of magnetic-field- 
induced pair-breaking effects in the superconducting component outside of the vortex cores. 

PACS numbers: 74.78.-w, 74.25.Ha, 74.25.N-, 74.25.nd 
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The vortex or Abrikosov state exists in type-II supercon- 
ductors subjected to magnetic fields between B cl and B c2 . 
In this state the field penetrates the superconductor in the 
form of quantized tubes of flux, or vortices. The supercon- 
ducting gap is zero inside the vortex cores and finite outside 
so that each vortex may be considered to have a core of nor- 
mal metal, surrounded by superconductor [1, 2]. Because 
vortex quantization renders the material an inhomogeneous 
system, it necessarily affects the electrodynamics of the su- 
perconductor. The microwave response of the vortex state 
has been extensively studied theoretically [3-5] and exper- 
imentally la- Hill . However the picture is still incomplete 
in the infrared region spanning the superconducting gap 

EH. 

In this Letter we address the infrared electrodynamics of 
the vortex state. We obtain the complex optical conductiv- 
ity of type-II superconductors and compare our results to 
calculations of a superconductor-normal metal mixture us- 
ing the two key models for the effective conductivity of an 
inhomogeneous system: that of Garnett IU4II (the so-called 
"Maxwell-Garnett theory" or MGT) and that of Brugge- 
man [15] (sometimes called the "effective-medium approx- 
imation" or EMA). We also compare our results to a the- 
ory of viscous motion of vortices driven by currents in the 
superconductor [3]. We find that only the MGT gives a 
good description of experiment, and then only when pair- 
breaking by the magnetic field Illa4l8ll is considered. That 
it does so is reasonable considering the topology of the 
vortex state: normal regions surrounded entirely by a con- 
nected superfluid. As pointed out some years ago [19], this 
is the topology of the MGT: the inclusions are embedded in 
a host medium and are correlated to stay apart. In contrast, 
the EMA allows percolation of the minority constituent at 
some critical concentration, something that does not hap- 
pen in the vortex state until the upper critical field, when 
the entire material is in the normal state. 

We studied type-II superconducting thin films of BCS 
superconductors Nbo.5Tio.5N and NbN, which are widely 
used in superconducting magnets |2fj], RF cavities Il2lll . 



and photodetectors [22]. The 10 nm Nbo.5Tio.5N film 
was grown on a quartz substrate in Ar and N 2 gas with a 
NbTi target, and the 70 nm NbN film grown on a MgO 
substrate in N 2 atmosphere using Nb, both by reactive 
magnetron sputtering 11231 UM. The substrates have neg- 
ligible absorption in the spectral range of interest (10- 
100 cm -1 ) for T < 20 K. We performed the experi- 
ment at Beamline U4IR of the National Synchrotron Light 
Source, Brookhaven National Laboratory. The beamline 
is equipped with a Bruker IFS 66v FT-IR spectrometer, 
modified to use synchrotron radiation and a superconduct- 
ing magnet for low-temperature magneto-spectroscopy. A 
composite silicon bolometer operating at T ~1.5K detects 
far-infrared radiation with high sensitivity. Both samples 
were cooled to 2 K (<C T c ) in zero field, and their trans- 
mission and reflection measured in magnetic fields from 
0-10 T, with the field direction normal to the sample sur- 
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FIG. 1. (color online) The vortex-state to normal-state trans- 
mittance ((a) and (c)) and reflectance ((b) and (d)) ratios for 
Nbo.5Tio.5N and NbN at selected perpendicular magnetic fields. 
The reflection data include a stray-light correction discussed in 
the text. 
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FIG. 2. (color online) The real (circles) and imaginary (squares) part of the vortex-state optical conductivity a v for Nbo.5Tio.5N and 
NbN, normalized to the normal-state conductivity a n . The lines are MGT fits, neglecting (solid lines) and including (dashed lines) the 
pair-breaking effects. 



face. To avoid systematic experimental errors, the samples 
were brought to the normal state by heating to 20 K in zero 
field without changing the position, and their transmission 
and reflection measured. These transmission and reflec- 
tion values served as a reference for calculating the real 
and imaginary parts of the optical conductivity relative to 
the normal state conductivity. An unavoidable reflection 
from a quartz vacuum window caused stray light to be in- 
cluded in reflection measurements [24]. This was corrected 
by measuring the stray light and subtracting. Spectra were 
collected at 4 cm -1 resolution, large enough to average out 
the interference fringes due to multiple internal reflections 
of the light in the substrate. The data for both samples at 
selected fields are shown in Fig. Q] As the magnetic field 
increases, the vortex-state transmission % and reflection 
1Z V continuously approach the normal-state values T n and 
lZ n . Nbo.5Tio.5N reverts to the normal state more quickly 
than NbN due to its lower B c2 . The samples were also 
cooled through T c at each field shown in Fig. Q] and the 
same measurements repeated. The results show no notice- 
able difference from the zero-field-cooled case. 

Measurements of both transmission and reflection en- 
able the direct extraction of the thin film optical conduc- 
tivity a = (j\ + iu%. We measured the absolute normal- 
state transmittance l~ n , using an open aperture as reference. 
The films have large electronic scattering rates 1/r that 
& n = °o/(l + iwr) ~ o"o in the far-infrared, yielding a 
frequency-independent T n . From the measured T n we ex- 
tract (T , which in turn determines lZ n . Using the method 
described in Ref. 11811 . we obtained o\ and <r 2 at various 
magnetic fields, shown respectively as circles and squares 
in Fig. |2] The extraction of a 2 is prone to errors at high 



frequency, because it involves a small difference of two 
relatively large numbers. <r 2 is further complicated by the 
residual fringes in the raw data above 40 cm -1 . The opti- 
cal conductivity data show two characteristic features: (i) 
With increasing magnetic field, the real part approaches the 
normal- state value (but does not exceed it whether below or 
above he gap frequency) while the imaginary part dimin- 
ishes, (ii) A weak minimum in the real part appears around 
the gap frequency at finite fields. 

We first compare these features with calculations based 
on two different effective medium theories for the average 
dielectric response of a mixture of two materials; the MGT 
and EMA 11141 fl5l l24l - l2ql . The nanometer dimension of 
a vortex core is much smaller than the far-infrared wave- 
length. This justifies treating the vortex-state superconduc- 
tor as an effective medium composed of a superconductor- 
normal metal mixture. The MGT and EMA make differ- 
ent assumptions, thus yielding different results. The MGT 
treats one constituent as the host and all others as embed- 
ded media, so is suitable for mixtures with isolated inclu- 
sions B19I1 - The effective response functions vary smoothly 
with the volume fraction / of the grains, and no percolation 
transition occurs for / < 1. In contrast, the EMA treats 
all constituents equivalently, appropriate for mixtures with 
connected grains. It predicts percolation at a certain mix- 
ture fraction f c . A vortex-state superconductor maintains 
its superconductivity even at a high volume fraction of vor- 
tices; percolation does not occur. Moreover, vortices cor- 
relate to stay apart because of the repulsive force between 
them. Therefore we expect the MGT to be a better descrip- 
tion of our experiment, though one previous infrared study 
argued otherwise IU3II . 
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FIG. 3. (color online) (a)-(c) Effective optical conductivity 
calculated from MGT (solid lines), EMA (dashed lines), and 
the Coffey-Clem model in the vortex-pinning regime at luq — 
300 cm -1 (dashed-dotted lines). / is the normal-volume fraction, 
(d) Effective optical conductivity calculated from the Coffey- 
Clem model at different depinning frequencies cjq for / = 0.4. 



Model calculation results in Fig. [3] support the above 
reasoning. When calculating the effective optical conduc- 
tivity, we treat the vortex cores as cylinders of dirty-limit 
Drude metal and the superconducting component using 
Mattis-Bardeen theory [27]. For the MGT the supercon- 
ducting component is taken as the host medium. Fig. Ha)- 
(c) compares the real and imaginary optical conductivity 
of the effective medium at various volume fractions / of 
the normal metal. The MGT results are consistent with the 
experimental results shown in Fig. |2] including the mini- 
mum in the real conductivity near the gap frequency and 
the monotonic growth toward a n for increasing field. The 
EMA results are strikingly different: the real conductiv- 
ity is greater than a n towards zero frequency, and changes 
non-mono tonic ally with field. Moreover, the 1/w behav- 
ior of the imaginary conductivity is completely suppressed 
in the EMA for / > 0.5, in contradiction to experiment 
(Fig. 12). 

The effective medium theories discussed above assume 
a static vortex lattice. To consider the contribution of 
vortex motion, we compare our data to the model pro- 
posed by Coffey and Clem [3], which calculates the surface 
impedance of type-II superconductors under the influence 
of vortex dynamics. They generalized the two-fluid model 
to couple self-consistently the supercurrent density with 
the vortex displacements. The central result of the model 
is a frequency-, field-, and temperature-dependent complex 



penetration depth X(uj,H,T) fHH]. Noting the relation 
between the complex penetration depth and the complex 
optical conductivity o"cc = */ '/"o^A 2 , and using the field- 
dependence of various quantities in Ref [28], we obtain the 
complex optical conductivity in the low temperature limit 
(T < T c ), 



ace _ (1 - b)a s /a n + b 
a n j3b(l -b)a s /a n + l" 



(1) 



b = B /B c2 is the reduced field. (3 = 1/(1 — iu /uj), 
where oj is the characteristic frequency that distinguishes 
the flux-pinning and flux-flow regimes, called the depin- 
ning frequency. a s is the optical conductivity of the super- 
conducting component. Eq. (Q]) reduces to correct limits as 
b — > and b — > 1. Since our superconducting samples 
were in the low temperature limit, we did not include the 
effect of thermal creep in deriving Eq. (Q3 and the expres- 
sion for p. 

Fig. (2d) compares the effective Coffey-Clem optical 
conductivity calculated from Eq. (Q]l at different uj . We 
assume b ~ / to allow comparison with the calculation re- 
sults of MGT and EMA, also shown in the figure. (The va- 
lidity of this approximation will be discussed below.) We 
note that using u = 300 cm -1 corresponds to the flux- 
pinning regime and uj = 3 cm -1 the flux-flow regime. In 
the flux-pinning regime, the real conductivity agrees with 
the data as well as with MGT and EMA above the gap, but 
almost flattens below the gap. Such behavior misses the 
minimum around the gap, found at different values of / in 
Fig. [3ta)-(c) and in the data in Fig. [2] (b,c,d,f,g,h). Flux 
flow brings a minimum to the real conductivity Re(<Tcc)> 
but it also substantially steepens the low-frequency part of 
Re(crcc) an d suppresses the high-frequency part. These 
features render the Coffey-Clem theory inconsistent with 
experiment. 

Having established that the MGT possesses the salient 
features of the experimental data, we now focus on a 
quantitative analysis. Previous studies of the vortex-state 
electrodynamics in superconductors assumed the super- 
conducting fraction to be unaffected by the magnetic field 
BLJJLlLJ]- Adopting this assumption, we use the Mattis- 
Bardeen theory 112711 to describe the superconducting frac- 
tion and the Drude model for the vortices. The volume 
fraction of the normal vortex cores / is varied to fit the 
real optical conductivity data at different fields. Once an 
optimal fit is found, we use / to calculate the imaginary 
part. The fits and calculations are shown as solid lines 
in Fig. |2j The zero-field case is trivial as we set / = 0. 
The minimum in the real conductivity around the gap fre- 
quency is captured by the fits, but appears shallower in the 
data than in the fits. This indicates less missing spectral 
weight in the data, meaning weaker superconductivity than 
assumed for the fits. The low-frequency imaginary conduc- 
tivity also corroborates this view. The calculated imaginary 
conductivity, especially that of NbN, is above the data as 



the field increases. Therefore, the analysis points to addi- 
tional weakening of superconductivity beyond the presence 
of the vortex cores. 

One source of such additional weakening could be 
magnetic-field-induced pair-breaking effects in the super- 
conducting fraction. This effect needs to be considered be- 
cause the magnetic field outside of the vortex cores does 
not diminish significantly until several penetration depths 
away from the core center [|29ll . The vortex spacing for a 
square lattice [30] is a\j = ^<& /B ~ 49 nm at 1 T and 
less at higher fields, where <J> is the magnetic flux quantum 
and B is the applied field. These distances are smaller than 
the penetration depth of our samples (a few hundred nm). 
Therefore the field penetrates all of the superconducting 
component for our measured field range. This field breaks 
the time -reversal symmetry of the electron pairing, result- 
ing in pair-breaking effects 11161, Il7ll . Such effects reduce 
the energy gap and smear out the gap-edge singularity in 
the quasiparticle density of states, consequently modifying 
the optical conductivity. These effects have been demon- 
strated in our studies of Nbn kTJq rN and NbN thin films in 
parallel magnetic fields IU8ll3lll . 

We re-evaluate the MGT analysis by including the pair- 
breaking effects. A complete model requires computing 
the field distribution in the superconducting fraction, and 
then the pair-breaking optical conductivity as a function of 
distance from the vortex center. For simplicity, we use the 
pair-breaking optical conductivity from our parallel-field 
study. Such an approximation is valid when the average 
field in the superconducting fraction is close to the applied 
field. In fact, for the field-normal-to-film geometry, tak- 
ing into account a demagnetization factor ~1, the average 
internal field almost equals the applied field IB2II . Taking 
the value of / from the previous fits, we calculate the effec- 
tive optical conductivity using MGT. The results are shown 
in FIG. Elb)-(d) and (f)-(h) as dashed lines. Inclusion of 
the pair-breaking effects significantly improves the quality 
of the fits, especially for NbN at high fields where such 
effects become significant. The calculated optical con- 
ductivity, both the real and imaginary parts, matches bet- 
ter to the data. The improvement is hardly noticeable for 
Nbo.5Tio.5N mainly because its B c2 at 2 K is about 11 T; 
thus the greater density of vortex cores dominates the aver- 
age response, leaving the pair-breaking effects difficult to 
see. 

The analysis presented above yields the normal-volume 
fraction /. Its value at different fields for both samples is 
plotted in Fig. |4] We previously assumed / = B / B c2 
when discussing the Coffey-Clem model. This assump- 
tion can be justified by noting that / = irrl/afj = 
TTr^B/^Q oc B, where we have assumed a square vor- 
tex lattice with lattice spacing an and vortex core radius 
r v . A linear fit yields B c2 = 10.9 T for Nbo.5Tio.5N 
and B c2 = 25.6 T for NbN. For comparison we deter- 
mined B c2 of both samples by four-probe resistivity mea- 
surements at Station SCM2 at the National High Mag- 



E 

o 
> 

75 
E 
o 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



1 1 1 
O Nb 05 Ti 05 N 


1 1 


" □ NbN 




Of 


-erl5n 


^^^ 1 1 1 


1 1 



2 4 6 8 

Magnetic field B [T] 



10 



FIG. 4. (color online) Field dependence of the normal-volume 
fraction / forNbo.5Tio.5N (circles) and NbN (squares). The solid 
lines are fits to / = B/B c2 . 



netic Field Laboratory [|24ll . We found B c2 ~ 11 T for 
Nbo.5Tio.5N and B c2 > 20 T for NbN. While the two 
methods give consistent result, we point out that the field 
dependence of / shown in Fig.|4]is sub-linear for both sam- 
ples. Such behavior is accounted for by the shrinking of 
the vortex cores with increasing field. Field-induced reduc- 
tion of the vortex core radius by more than a factor of two 
has been observed in conventional s-wave superconductors 
[30], explained in terms of either the "vortex lattice squeez- 
ing effect" or the inter- vortex transfer of bounded quasipar- 
ticles in the vortex cores. Although the latter effect was ar- 
gued to be significantly weakened by scattering [33], it has 
been observed in dirty Pb thin films B4II . When the vortex 
core radius r v (B) decreases with field, the normal-volume 
fraction / = irrl(B)B/& could become sub-linear in 
field. 

In conclusion, we measured the far-infrared transmission 
and reflection of Nbo.5Tio.5N and NbN thin films in per- 
pendicular magnetic fields, and observed superconductiv- 
ity weakening. The extracted optical conductivity of the 
vortex state is consistent with the Maxwell-Garnett theory 
for normal metal inclusions in a superconducting host. De- 
tailed analysis suggests that in the vortex state the magnetic 
field weakens superconductivity via at least two ways. One 
is the increasing normal metal component as the vortex 
density increases with field. The other is the weakening 
of the BCS state in the remaining superconducting com- 
ponent due to pair-breaking effects. We did not need to 
include the flux-flow-induced dissipation in our analysis. 
Effects of flux flow are expected to be important in cleaner 
materials with less pinning. 
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Supplemental Material 

I. Sample Parameters 

The sample parameters relevant to our study are listed 
in Table |Sl] The film thickness d was estimated from the 
film growth conditions. The zero-temperature gap A was 
extracted from the zero-field optical conductivity, shown in 
Fig. 2 of the main text. The critical temperature T c was de- 
termined from temperature-dependent resistivity measure- 
ments. The normal-state sheet resistance Rq was deter- 
mined from the normal-state transmittance, 

_ An 

'" (n + l + V^n) 2 ' { } 

where n is the substrate refractive index and Z is the vac- 
uum impedance. We found the penetration depth A, co- 
herence length £, and the substrate refractive index n from 
literature yj-|3|]- 

To determine the upper critical field B c2 of both sam- 
ples, we performed four-probe resistivity measurements at 
Station SCM2 at the National High Magnetic Field Lab- 
oratory. A 2 mmx2 mm piece was cut from each sam- 
ple for such measurements because of the requirements of 
the experimental set-up. Electrical contacts were made by 
gluing gold wires to the thin films with silver paint. The 
two samples were mounted simultaneously on a rotating 
probe, which allows accurate alignment of the sample sur- 
face with respect to the magnetic field orientation. Resis- 
tivity was measured with the field normal to the film sur- 
face. Well below the critical temperature, the resistance 
was measured by scanning the magnetic field between T 
and 16 T at selected temperatures. Close to the critical tem- 
perature, the resistance was measured by sweeping tem- 
perature between 5 K and 20 K at selected low field val- 
ues. The data are shown in Fig. [STj Because the transi- 
tion is broad, we define B c2 as the field at which the re- 
sistance drops to 1/2 of its 20 K value. We determined 
B c2 « 11 T at 2 K for Nb . 5 Ti . 5 N. B c2 at 2 K for NbN 
is beyond the measurement range of the magnet. Based 
on the temperature-dependent trend shown in Fig. |Sl] we 
estimate that it should be greater than 20 T at 2 K. 

II. Infrared Measurement Procedures 

For transmission and reflection measurements, the sam- 
ples were loaded into the superconducting magnet at the 





TABLE SI. 


Sample parameters 


Sample 


d 
(nm) ( 


Ao 
—i 
cm 


> (K) 


Ra 


A £ n 
(nm) (nm) 


Nb . 5 Ti . 5 N 
NbN 


10 
70 


12.8 
17.9 


10.2 
12.8 


117 200-400 3.8-5.0 2.12 
48 180-500 4.0-7.0 2.90 



bottom of a probe in the Faraday configuration (external 
magnetic fields perpendicular to the film surface). We 
started with transmission measurement. In the experiment 
layout shown in Fig.[S2] the plane mirror (A) for collecting 
the reflected light was taken out from the four-way cross- 
shaped flange, so that the incident beam could completely 
pass through. The transmitted signal was detected by the 
bolometer placed at Position PI. To apply the field perpen- 
dicular to the film, we first roughly loaded the samples in 
the approximate angle range, and then carefully rotated the 
probe to maximize the transmission signal. The position 
was carefully marked for later reference. The error of the 
alignment is expected to be within 5°. 

The reflection measurement was done immediately af- 
ter the transmission measurement. The sample was kept 
at the same position, with the orientation of the sample 
probe unchanged so that the magnetic field was perpen- 
dicular to the film surface. Plane mirror A was installed to 
direct the reflected signal from the sample to the bolome- 
ter, which had been moved from Position PI to Position 
P2 downstream of the plane mirror. There are two sets of 
quartz windows, one set on the exterior of the magnet, the 
other set on the Variable Temperature Insert (VTI) in which 
the sample was mounted. The stray light was mainly from 
the entrance quartz window on the VTI. This stray light 
was measured by rotating the sample probe 45° from the 
marked position mentioned in the previous paragraph, and 
subtracted from all the single beam reflection spectra of the 
samples. To minimize the stray light, the plane mirror (A) 
and the paraboloid (B) were manipulated to minimize the 
ratio of the stray signal and the signal from the sample at 
the marked position (true signal plus stray signal). 



III. Effective Medium Theories 

Here we derive the effective optical conductivity of a 
vortex-state superconductor from the Maxwell-Garnett the- 
ory and the Bruggeman effective medium approximation. 
The basic formulation of effective medium theories can be 
found in literature, e.g. Q. 

Consider an inhomogeneous medium made of two com- 
ponents, grain a with volume fraction / embedded in the 
surrounding b with volume fraction 1 — /. The electro- 
magnetic waves traversing the medium can be treated as 
a spatially-averaged field with its electric component ex- 
pressed as 



(E) = /E a + (1 - /)E 6 



(S2) 



The response function D is similarly spatially averaged, 

(D) = fe a B a + (1 - f)e b B b . (S3) 

The effective dielectric function is then 
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FIG. SI. Four-probe resistance data in perpendicular fields for Nbo.5Tio.5N (first row) and NbN (second row). Left column: resistance 
vs perpendicular field at different temperatures. Middle column: resistance vs temperature at different perpendicular fields. Right 
column: circles are B c2 extracted from data. In the resistance plots dashed lines indicate the resistance at 20 K and half of that value. 
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FIG. S2. Optical layout for transmission and reflection measurements in a perpendicular magnetic field. The red rays are the incident 
beam, the green rays the transmitted beam, and the blue rays the reflected beam. A plane mirror is installed in the four-way cross-shaped 
flange to direct beams in the reflection measurement, and is taken out to allow the full beam to pass in the transmission measurement. 
The bolometer is placed at Position PI for transmission and moved to Position P2 for reflection measurement. 



1. Maxwell-Garnett theory 

The Maxwell-Garnet theory considers the electric field 
E a to be the local field acting on the grain a. It consists 
of the external field and the field caused by the polarized 
charges on the surface of an artificial cavity in which the 
grain resides. In considering this local field, the theory as- 
sumes that the separation between grains are large enough 
so that individual grains scatter light independently. Fur- 
thermore, the field in the surrounding medium b is assumed 
to be unaffected by the presence of the grains 0.0]. Based 
on these assumptions, the Maxwell-Garnett theory gives 
the effective dielectric function for oriented ellipsoid grains 
as 



£mgt — e& + £fc 



/(e a - e&) 



9(1 - /)(e a -e b ) + e b 



(S5) 



where / is the volume fraction of grain a, and g is the 
depolarization factor determined by the shape of the el- 
lipsoid. To convert to the optical conductivity, note e = 
1 + Airier /u, and a a = a n and o b = a s . Vortices are gen- 
erally regarded as cylindrical tubes with a normal core of 
radius of the coherence length £, each carrying a quantum 
of magnetic flux <J> . Therefore we set g = 1/2. 



2. Bruggeman effective medium approximation 

Bruggeman proposed a method to address the issue that 
the grains and host material in Maxwell-Garnett theory are 
treated asymmetrically. Because of the presence of grains 
with different properties from the surrounding medium, the 
electric field in the region around the grains is modified, 
and the electric flux deviates from that when such grains are 
absent. Bruggeman argued that the average flux deviation 
for the whole medium should vanish J3]. He then suggested 
that an adequate choice of a self-consistent local field can 
satisfy this condition. This leads to the consideration of 
an effective medium in which all inclusions are treated on 
an equal basis, and the average flux deviation is zero [|6J]. 
Equivalently, the effective dielectric function can be calcu- 
lated using Eqs. ([S2]), ([S3]l, and ([S4]). The difference from 
the Maxwell-Garnett theory is that there is no host medium 
in Bruggeman effective medium approximation. Both E a 
and E h are calculated assuming the constituents a and b 
are grains immersed at their related volume fractions in an 
effective medium characterized by an effective dielectric 
function. 



For an inhomogeneous medium made of two compo- 
nents, grain a of volume fraction / and grain b of volume 
fraction 1 — /, Bruggeman effective medium approxima- 
tion gives the effective dielectric function for oriented el- 
lipsoid grains as the solution of the following equation, 



/ 



e EMA 



ge a + (1 - sOcema 



+ (1-/) 



e fc — ^EMA 



ge b + (1 - g)e E MA 



= 0. 
(S6) 



Treating the vortices as cylindrical particles, we set g 
1/2. The solution is 



1 



£ema 



(2/-l)(e a -e 6 ) 



(2/-l) 2 (e a -e b ) 2 + 4e Q e b 



(S7) 



Only this one out of the two solutions of the quadratic equa- 
tion JS6I 1 is chosen, because it reduces to the correct values 
in the limit / — > and / — > 1. The optical conductivity is 
obtained using the relation e = 1 + Airier /uj, and a a = a n 
and a b = a s . 
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